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In this work we investigate the vahdity limits of the modulational approximation as a method to 
describe the nonlinear interaction of conservative wave fields. We focus on a nonlinear Klein-Gordon 
equation and suggest that the breakdown of the approximation is accompanied by a transition to 
regimes of spatiotemporal chaos. 
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I. INTRODUCTION 

Modulational instability of high-frequency nonlinear 
waves is a common process in a variety of circumstances 
involving wave propagation in continuous systems. Mod- 
ulational processes can be seen to occur in a wide range 
of physical situations, from nonlinear waves in plasmas 
to nonlinear electromagnetic waves propagating in op- 
tical fibers What usually happens in all those cases 
is that due to generic nonlinear interactions, the ampli- 
tude of a high-frequency carrier develops slow modula- 
tions in space and time. If the modulations are indeed 
much slower than the high-frequencies involved, one can 
obtain simplified equations describing the dynamics of 
the slowly varying amplitudes solely, the amplitude equa- 
tions In the present analysis we consider systems 
that become integrable in this modulational limit, a fea- 
ture often displayed. Should this be indeed the case, no 
spatiotemporal chaos would be observed there. The ba- 
sic interest then would be to see what happens when the 
approximations leading to modulational approximations 
cease to be satisfied. The paper is organized as follows: in 
§2 we introduce our model equation and discuss how and 
when it can be approximated by appropriate amplitude 
equations; in §3 we investigate the modulational process 
from the point of view of nonlinear dynamics; in §4 we 
perform full spatiotemporal simulations and compare the 
results with those obtained in §3; and in §5 we conclude 
the work. 



II. MODEL EQUATION, MODULATIONAL 
APPROXIMATIONS, AND AMPLITUDE 
EQUATIONS 

In the present paper we focus attention on a nonlinear 
variant of the Klein-Gordon equation (NKGE) to inves- 
tigate the breakdown of modulational approximations in 



the context of nonlinear wave fields, 
here reads 



The NKGE used 



d'iA{x,t)^dlA(x,t) -f 
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{dt = d/dt, dx = d/dx) where we write the generalized 
nonlinear potential as 
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Lo playing the role of a linear frequency which will set the 
fast time scale. The remaining coefficients on the right- 
hand-side were chosen to allow for modulational instabil- 
ity and saturation; we shall see that while the negative 
sign of the second fulfills the condition for modulational 
instability, the positive sign of the third provides satu- 
ration. The choice of their numerical values is arbitrary, 
but our results are nevertheless generic. The NKGE is 
known to describe wave propagation in nonlinear media 
and the idea here is to see how the dynamics changes as 
a function of the parameters of the theory: wave ampli- 
tude, and time and length scales. 

Let us first derive the conditions for slow modulations. 
We start by supposing that the field A{x, t) be expressed 
in the form 



A{x,t) = A(x,f)e'"* -f complex conjugate. 



(3) 



Then, if one assumes slow modulations and discards 
terms like d^A and the highest-order power of the po- 
tential one obtains 

2iLudtA{x,t) - dlA{x,t) - 3\A{x,t)fA{x,t) = 0, (4) 

which is, apart from some rescalings, the Nonlinear 
Schrodinger Equation - we shall refer to it as NLSE here 
- an integrable equation. 
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This is no novelty; it is known that the modulational 
approximation is obtainable when there is a great dispar- 
ity between the time scales of the high-frequency uj and 
the modulational frequency, we call it SI, such that terms 
of order Sl^A, when compared to uP'A, can be dropped 
from the governing equation. The magnitude of the mod- 
ulational frequency can be estimated as follows. Con- 
sider Eq. and suppose a democratic balance among 
the magnitude of its various terms; LodtA ~ ff^A ~ . 
Then one obtains for dt —> fi, 



(5) 



It is thus clear that the modulational approximation is 
valid only when A <^ ld^ since this condition slows down 
the modulational process causing f2 w. The next ques- 
tion would be on what is to be expected when the mod- 
ulational approach ceases to be valid. Before proceeding 
along this line, let us mention that a stability analysis 
can be performed on Eq. (^). One perturbs an homo- 
geneous self-sustained state with small fluctuations of a 
given wavevector fc > (we choose fc > here, but the 
theory is invariant when fc — > — fc) and after some alge- 
bra one concludes that: (i) the field in the homogeneous 
state, let us call this field Ah, is given by 



Our first task is to examine how the regular and well 
known modulational instability analyzed in the previ- 
ous section comes directly from Eq. (|^). To do that, 
let us write a truncated solution as the sum of an ho- 
mogeneous term plus fluctuations with wavevector fc, 
A = Ah{t) + Ai(i)(e*'=^ -I- e-^*^^). Ah and Ax real. The 
truncation, that discards higher harmonics, is legitimate 
within linear regimes, but fortunately we shall see that it 
is not as restrictive as it might appear even in nonlinear 
regimes. The general idea favouring truncation here is 
that in the reasonable situation where modes with the 
fastest growth rates are more strongly excited, relation 
indicates that second harmonics are already outside 
the instability band since 2kmax — V^ktr > ktr- Under 
these circumstance one would be led to think that the 
most important modes would be the homogeneous and 
those at the fundamental spatial harmonic. We shall ac- 
tually see that the truncation provides a nice and repre- 
sentative approach to the case of regular regimes. 

Note that we are already considering the amplitudes of 
the exponential functions as equal. This results from a 
simplified symmetrical choice of initial conditions and is 
totally consistent with the real character of Eq. (|^). Af- 
ter some lengthy algebra, one finds out that the coupled 
nonlinear dynamics of the fields Ah and Ai is governed 
by the Hamiltonian 



Ah = flc 



(6) 



where is an arbitrary amplitude parameter and where 
the exponential term should be seen as providing a small 
nonlinear correction to the linear frequency uj; and (ii) 
the perturbation is unstable when 



k <ktr = \/6a 
with maximum growth rate at 

V2 



(7) 



(8) 



When unstable, the homogeneous state typically evolves 
towards a state populated by regular structures which 
can be formed precisely because the underlying governing 
NLSE, Eq. (|), is of the integrable type as mentioned 
earlier. 



III. BEYOND THE MODULATIONAL 
APPROXIMATION 



To advance the analysis beyond the modulational 
regimes we start from the basic equation, Eq. (]^), but 
do not use approximation d^A <C oj'^A leading to Eq. (Q) 
and eventually to condition A ^ lu. The idea is precisely 
to examine what happens as the ratio A/oj of Eq. (1q) 
grows from values much smaller, up to values compara- 
ble to the unit. 
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where = t^^+fc^, 9h = Ah/\/2, qi = Ai, and where the 
p's denote the two momenta conjugate to the respective 
g-coordinates. 

The Hamiltonian (^) can be informative. As a first in- 
stance it can be used to determine the stability properties 
of the homogeneous pump, as mentioned before. To see 
this, assume that in average qh ^ qi and solve the dy- 
namics perturbatively. In zeroth order, one would have 
the following Hamiltonian governing the (ph^Qh) dynam- 
ics: 



^° " 2 + 2 2 + 3 ■ 



(10) 



With help of action-angle variables (p, O) for the zeroth- 
order part and conventional perturbative techniques the 
solution reads 



qh = \ — cos(a;i- -z^t)^ 



(11) 



if yO, the amphtude parameter, is not too large. Note 
that the oscillatory frequency undergoes a small nonlin- 
ear correction which is determined by the quartic term 
in qh of the Hamiltonian ho. 

Next we consider the driven Hamiltonian controlling 
the dynamics of the canonical pair (pi, q\) 



2 



hi 
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where we recall that the pair (pi , qi ) describes the inho- 
mogeneity of the system. One again introduces action- 
angle variables (/, 6) to rewrite the linear Hamiltonian 
(|l|) in the for m 



hi = 12--cos2 6lcos^M- ^-^t), 
from which we obtain the resonant form 



hi^r 



1 
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^cos(20), 



(13) 
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with (/) ~ — [uj — 3p/2uj'^)t, and where use is made of 
the approximation x ^ uj -\- fc^/(2a;) valid when k uj. 
Now consider the stability of a small perturbation / 
0. For this initial condition hi_r 0. Since hi^r is a 
constant of motion, solutions for arbitrarily large values 
of /, what would indicate instability, are possible only 
when |cos(2(/i)| < 1. From the resonant Hamiltonian 
(p!^), this demands 



(15) 



and also indicates that maximum growth rate for / occurs 
at kmax = y^ph^- Comparisons of temporal dependence 
of Eqs. (H), (H) and (|ll| ) shows that pjio = a^, and that 
conditions ^ and (|15D are therefore one and the same 
as they should be. In other words, starting from the full 
nonlinear wave equation, we recover the typical results 
naturally yielded by the NLSE. 

But the Hamiltonian (^ gives further information be- 
cause apart the truncation, it is not an adiabatic approx- 
imation like Eq. (^; it can therefore tell us whether the 
reduced dynamics, if unstable, is of the regular or chaotic 
type. The interest on this issue comes from the fact that 
the reduced dynamics usually helps to determine the spa- 
tiotcmporal patterns of the full system: while regular 
reduced dynamics is associated with regular structures 
- frequently a collection of solitons or soliton-like struc- 
tures, chaotic reduced dynamics is associated with spa- 
tiotemporal chaos. The correlation has its roots on the 
so called stochastic pump model The model states 
that intense chaos in a low-dimensional subsystem of a 
multidimensional environment can make the subsystem 
act like a thermal source, irreversibly delivering energy 
into others degrees-of-freedom. In the limit of a predom- 
inantly regular dynamics undergoing in the reduced sys- 
tem, irreversibility is greatly reduced and energy tends 
to remain confined within the subsystem. On the other 
hand, in the limit of deeply chaotic dynamics with no pe- 
riodicity at all, energy flow out of the subsystem is fast 
and there is not even much sense in defining the subsys- 
tem as an approximately isolated entity. The interme- 
diary cases are those more amenable to a description in 



terms of the stochastic pump (|,^. We point out, how- 
ever, that in all cases, even in the chaotic one, analysis of 
the reduced subsystem serves as an orientation on what 
to be expected in the full spatiotemporal dynamics. In 
our case, the appropriate subsystem is precisely the one 
we have been using. This is so because it is the smallest 
subsystem comprising the most important ingredients of 
the dynamics: the homogeneous and the only linearly 
unstable modes. 

We now examine these points in some detail. Let us 
consider lo — 0.1. When A^^o <C w, A^^o = Ah{t — 0), 
the modulational approximation can be used to deter- 
mine the instability range, k < ktr- A surface of sec- 
tion plot based on the two-degrees-of-freedom Hamilto- 
nian (^ produces Fig. (Q), where we record the values 
of (qi,pi) whenever = with dph/dt > 0, and where 
we take qi — O.Olqh ^ qh and ph = pi = to determine 
the total unique energy of the various initial conditions 
we launch in the simulations. Note that the particular 
"seed" initial condition introduced above is included in 
the ensemble of initial conditions launched and represents 
small perturbations to an homogeneous background. 
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FIG. 1. Low-dimensional Poincare plots on the projected 
phase-space u — 0.1 in all cases, and k — kmax in 

panels (b)-(d). Initial conditions discussed in the text. 

In Fig. (|l]a) we take Ah^o/i^J = 0.5 and k = 2kmax > 
ktr- One lies outside the instability range and the fig- 
ure reveals that the origin, which represents the purely 
homogeneous state qi = Ai = 0, is indeed stable. The re- 
maining panels are all made for k = kmax and increasing 
values of the ratio Ah.o/t^- One sees that for such values 
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of k and A^^o not only the origin is rendered unstable, 
as it also becomes progressively surrounded by chaotic 
activity. The inner chaotic trajectory issuing from the 
origin - i.e., the trajectory representing the modulational 
instability - is encircled by invariant curves, but the last 
panel already shows that some external chaos, here repre- 
sented by scattered points around the invariant curves, is 
also present. At some amplitude Ah^o = Acriucai slightly 
larger than the one used in panel (d), invariant curves 
are completely destroyed. Above the critical field, inner 
orbits are no longer restricted to move within confined re- 
gions of phase-space - these orbits are in fact engulfed by 
the external chaos seen in panel (d) . External chaos here 
is a result of the hyperbolic point positioned at the local 
maximum of the generalized potential ^{A), at Ai ^ 0.1. 
One gross determination of the critical field based on sim- 
ple numerical observation yields AcriUcai / 1-^ ^ 1/1.62, al- 
though the corresponding transition from order to disor- 
der in the full simulations may not be so sharply defined. 
Yet, one may expect that in both low-dimensional and 
full simulations, the transition should occur at compara- 
ble amplitudes. 



IV. FULL SPATIOTEMPORAL SIMULATIONS 



To make the appropriate comparisons, we now look 
into the full simulation of Eq. (|l|). Full simulations are 
made through a discretization of the spatial domain via a 
finite difference method. The dynamics is resolved tem- 
porally by means of a symplctic integrator, as was the 
purely temporal Hamiltonian. The results are quite ro- 
bust and energy is conserved up to 10~^ parts in one. In 
Fig. (H) we display the spacetime history of the quantity 



Q — \J uj'^A{x,tY + A{x,tY - we plot this quantity be- 
cause it becomes a constant of motion in the limit where 
we discard nonlincarities and inhomogeneities. Initial 
conditions are the same as the seed initial condition used 
in Fig. dl), but the control parameters differ. In the case 
of Fig. (ga) where we choose Ah.o/^^ — 0.1 so as to safely 
satisfy A-^^o <C Acriucai and Ah^o ^ it is seen that the 
spatiotemporal dynamics is very regular. The homoge- 
neous state is unstable, but only periodic spatiotemporal 
spikes can be devised. This is the regular spatiotemporal 
dynamics so typical of the integrable NLSE and this kind 
of dynamics agrees very well with the low-dimensional 
predictions of the reduced Hamiltonian (^. The fact 
that only one single structure can be found along the 
spatial axis at any given time, indicates that the dynam- 
ics is singly periodic (period- 1) along this axis and thus 
basically understood in terms of the reduced number of 
active modes (homogeneous plus fundamental harmon- 
ics) of Hamiltonian (|^). We now move into the vicinity 
of the critical amplitude AcriUcai discussed earlier. Un- 
der such conditions, one may expect to see the effects 
of spatiotemporal chaos. The value A/i o/w = 1/1.625 is 
chosen in Fig. where full simulations indeed display 



a highly disorganized state after a short regular transient. 
Regularly interspersed spikes can no longer be seen and 
spatial and temporal periodicities are lost, which charac- 
terizes spatiotemporal chaos. In this regime many modes 
become active (we shall return to this point later) and 
the reduced Hamiltonian (^ fails to provide an accurate 
description of the full dynamics. However it still pro- 
vides a good estimate on the point of transition. A more 
throughful examination of the transition in the full simu- 
lations suggests that the critical field there is a bit smaller 
- a value close to a;/1.95; for smaller values we have not 
observed noticeable signals of spatiotemporal chaos even 
for much longer runs than those presented here. The 
much larger oscillations executed by Q in chaotic cases 
(see the legend of Fig. (||)) is a direct result of the de- 
struction of the invariant curves seen in Fig. (0). In 
the absence of invariant curves, initial conditions are no 
longer restricted to move near the origin. 




spatial grid 



(b) 



spatial grid 

FIG. 2. Full spatio-temporal simulations for lj = 0.1 and 
k = fcmaj;. Ah,ol'^ = 0.1 in (a) and Ah,ol^ = 1/1.625 
in (b). The intensity plots are made for the quantity 

Q = \/ ixP- A(x ^tY' 4- ^(x, t)2; lighter shades are associated 
with larger values of Q, which varies within the range 
< Q < 0.0014 in (a) and < Q < 0.8 in (b). t = ujijl-K 
here and in all the remaining simulations 

It is thus seen that there are limits to an integrable 
modulational description of the dynamics of a wave field. 
The limits are essentially set by the parameter A/j q/oj. 
If it is much smaller than the unit, the modulational 
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description is valid and one can expect to see a collec- 
tion of spatiotemporal periodic structures being formed 
as asymptotic states of the dynamics. On the other hand, 
as the parameter approaches the unit, nonintegrable fea- 
tures are likely to be seen. In particular, regularity does 
not survive for very long, and fluctuations with various 
length scales appear in the system. This is a regime of 
spatiotemporal chaos which fundamentally involves the 
presence of nonlinear resonances between the frequency 
of the carrier, uj, and the intrinsic nonlinear modulational 
frequency, fl. 

Due to the presence of chaos, one is suggested that the 
transition involves an irreversible energy flow out of the 
reduced subsystem. If one computes the average number 
of active modes 



< >E 



EM' 



(16) 



where the amplitudes are defined in the form 

^„ = ^A(x„^)e"'="^ (17) 

i 

one obtains the plots shown in Fig. 



A- 




V<N=> , 



FIG. 3. Average number of active modes as a function 
of time. Parameters respectively equal to those of the 
previous figure with exception of panel (d) where we take 
Ah,o/Lu= 1/1.95. 



In Eq. (|17|), "j" is the discretization index. In Fig. 
(^) we use the same conditions as in Fig. (^). This 



panel shows that in rough terms, energy keeps period- 
ically migrating between the homogeneous mode (when 
\/<N^> ^ 0) and the fundamental harmonic (when 
V< iV2 > - 1). Conditions of Fig. (|b) are the same 
as those of Fig. (|b); one sees that as the ratio Ah.o/^ 
grows, periodicity is lost, and that energy flow out of the 
initial subsystem into other modes becomes clearly irre- 
versible. In Fig. ^p) we use the same previous conditions 
with exception of A^^o which we take A^^o — uj/1.95. 
This slightly smaller, but not too small, value of the ini- 
tial amplitude allows to observe the slow diffusive transit 
of energy during the initial stages of the corresponding 
simulations. During this stage one can actually look at 
the subsystem as an energy source adiabatically deliver- 
ing energy into other modes - the concept of the stochas- 
tic pump applies more appropriately in those situations. 



V. FINAL CONCLUSIONS 

To summarize, in this paper we have studied the break- 
down of modulational approximations in nonlinear wave 
interactions. We have analyzed a nonlinear Klein-Gordon 
equation to draw the following conclusions. Adiabatic 
or modulational approximations are accurate while the 
high-frequency of the carrier wave keeps much larger than 
the modulational frequency. Under these circumstances 
the full spatiotemporal patterns are regular as is the dy- 
namics in the reduced subsystem where the system en- 
ergy is initially injected. There is no net flow of energy 
out of the reduced subsystem into the remaining modes. 

On the other hand, when both frequencies become of 
the same order of magnitude, the reduced subsystem un- 
dergoes a transition to chaos. Correspondingly, the spa- 
tiotemporal patterns of the full system become highly 
disordered and energy spreads out over many modes. 
The correlation between the low-dimensional and high- 
dimensional spatiotemporal chaos has its roots on the 
stochastic pump model |3|. According to the model, a 
low-dimensional chaotic subsystem may act like a ther- 
mal source, delivering energy in a irreversible fashion to 
other degrees-of-freedom of the entire system. Spectral 
simulations performed here indicates that this seems to 
be the case with the present setting. 

The transition to chaos involves a noticeable increase 
in terms of wave amplitude. This takes place when the 
invariant curves of Fig. (^) are destroyed, allowing for 
the merge of the external and internal chaotic bands into 
one extended chaotic sea. This merging of chaotic bands 
is actually a result of reconnections involving the mani- 
folds of the unstable fixed point at the origin, and other 
hyperbolic points associated with the curvature of the 
generalized potential $ ||]. When both manifolds recon- 
nect, inner trajectories issuing from the origin start to 
execute the large and irregular oscillations seen in Fig. 
(^). More detailed studied of the process is under cur- 
rent analysis. 
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